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Abstract 

We construct the wave functions for the Moore- Read v = 5/2 quantum Hall state on a torus in 
the presence of two quasiholes. These explicit wave functions allow us to compute the monodromy 
matrix that describes the effect of quasihole motion on the space of degenerate ground states. 
The result agrees with the recent discussion by Oshikawa et al. (Annals of Physics 322 1477). 
Our calculation provides a conformal field theory explanation of why certain transitions between 
ground states are forbidden. It is because taking a quasihole around a generator of the torus can 
change the fusion channel of the two quasiholes, and this requires a change of parity of the electron 
number in some of the ground states. 

PACS numbers: 73.43.Fj, 73.43. Jn 
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I. INTRODUCTION 

The degeneracy, first noted in numerical work by Yoshioka et al. and W. P. Su of 
the fractional quantum Hall ground state on torus is an essential feature of such states. 
In its absence, Laughlin's general gauge argument would require the Hall conductance to 
be an integer. Explicit wave functions for Laughlin states on torus were written down by 
Haldane and Rezayi [2J, who identified the degeneracy as arising from translations of the 
center of mass of the incompressible Hall fluid. Later, Wen and Niu realized [3( that the 
degeneracy persists even in the absence of translation symmetry, and that the degree of 
degeneracy was sensitive to the global topology of the space in which the Hall fluid resides. 
This observation lead them to introduce the notion of topological order as a characterization 
of the strongly-correlated ground state. Moore and Read 4| revealed the more general 
structure of topological order by pointing out its connection to rational conformal field 
theory. Since topological order is distinguished by the quantum numbers of ground state and 
excitations, a change of topological order requires a quantum phase transition corresponding 
to a substantial rearrangement of the many-body Hilbert space. Such phase transitions are 
not usually associated with symmetry breaking, and the resulting topological phases posses 
no conventional order parameters. They therefore fall outside the conventional Landau 
theory of phase transitions. 

It is possible for a quantum fluid to have topologically degenerate states even when 
it lives on the plane. The degree of degeneracy then depends on the number and type 
of vortex-like defects in the fluid. Braiding these defects produces transitions between the 
topologically protected degenerate states, and it has been suggested that such manipulations 
may be exploited for quantum computation [5]. One such topological phase, the Moore- Read 
Pfaffian state, is the likely candidate for the observed v = 5/2 quantum Hall state \&\. The 
Pfaffian wave function is easy to write down on the plane. When the fluid is placed on 
a torus, however, the wave function becomes more complicated. Counting the number 
of degenerate states is harder than counting the number states for a Laughlin fractional 
quantum Hall state, and the answer depends on whether the number of electrons in the 
system is odd or even. Nonetheless, the ground state wave functions have been constructed 
I fl a - although some issues still remai , U3h, g geaera, pri nc ipl e s that do aot require 
knowledge of the explicit Pfaffian wave functions, and by building on what is understood for 
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Abelian quasiparticles [9|, Oshikawa et al. [10|| have identified the topological quasi-particle 
operations that transform one Moore-Read ground state to another. 

In this paper we investigate how the topological operations of Oshikawa et a/.'s affect 
the many-electron wave functions. To do this, we first construct the Moore-Read Pfaffian 
states on a torus in the presence of two quasihole excitations. We then explicitly exhibit 
their transformation properties under topological operations on the quasihole positions. In 
this process, we provide a conformal field theory explanation of why the the monodromy 
matrix acting on the space of the degenerate ground states should be block-diagonal with 
respect to the parity of the spin structure. It is because the spin structure dictates both the 
fusion channel of two quasiholes and the parity of the number of electrons. 

II. LAUGHLIN STATES 

Before tackling the Pfaffian state, it helps to recall the properties of the simpler Laughlin 
states on a torus {2], Q, 11]. A torus can be regarded as a rectangle with sides L x and 
L y periodic boundary conditions. In the Landau gauge, with A = —By it, the periodic 
boundary conditions on the wave function are twisted by a gauge transformation that is 
necessary to prepare a particle that leaves the top edge of the rectangle for its reappearance 
at the bottom edge. For a single electron we have 

ip(x + L x ,y) = ?p(x,y) ij){x, y + L y ) = exp(—iL y x/l 2 )ijj(x, y), (1) 

where / = \JhjeB is the magnetic length of the system. This twisted boundary condition 
must continue to hold, particle by particle, for the many-particle wave function. In our 
Landau gauge, the lowest Landau level many-particle wave function takes the form 

(ri, r 2 , . . . , r Ne ) = exp(- ^ y 2 /2l 2 )f{z u z 2 ,...,z Ne ), (2) 

i 

where Z{ = (x, + iyi)/L x and / is a holomorphic function in each of the Zi. Combining 
Eq. (JTJ and Eq. (j2J) shows that / must satisfy the following quasi-periodicity conditions: 

f(z!,...,Zi+ 1,...) = f(zi,...,Zi,...), 

f(z 1 , . . . , Zi + r, . . .) = exp[-mN s (2 Zi + r)]./(zi, . . . , z h . . .), (3) 

Here r = iL y /L x and N s = L x L y j2 r nl 2 is the number of flux quanta passing through the 
torus. (These conditions on / preserve their form for a torus which is a periodic parallel- 



ogram, rather than a rectangle. In this case r is no longer purely imaginary, but should 
retain a positive imaginary part) 

By integrating d/dzi[\ia f(zi, . . . , ZjvJ] around the boundaries of the rectangle, and using 
Eq. ([3]) to combine the contributions of the opposite sides, we see that the number of zeros of 
/, considered as a function of Zi, is N s . The precise locations of these N s zeros will depend 



on the positions of the the other Zj, but they are subject to a non-trvial constraint 11]. 
Suppose that g(z) is meromophic and doubly-periodic: 

g{z) = g{z + l) = g{z + r). (4) 

By evaluating the integral 

< 5 > 

around the edges of the period parallelogram + — > 0, and again combining 

the contributions from opposite sides, we find that 

I = m + nr (6) 

where m and n are integers. On the other hand, if the poles of g(z) in the period parallel- 
ogram are at z = b,i and the zeros at z — then / = ^\ — hi. Taken together with 
Eq. (161), this means that the sum of zeros minus the sum of poles vanishes modulo periods. 

n 

This is Abel's theorem for the torus [12]. Although f(zi, . . . , Zjsrj is it not doubly periodic, 
the ratio 

g( z ) = f( z ,,z 2 ... } z Ne )/f(z } z' 2 ... } z' Ne ) (7) 

is a doubly-periodic meromorphic function whose zeros are at the zeros of f(z,Z2-.-, -Zjv e ) 
and whose poles are at the zeros of f(z,z' 2 . . . , z' N ) . Consequently the sum of N s zeros of 
/(/Zi, -22, ... , znJ considered as a function of Z\ is independent of other electron coordinates. 
The same is true of the zeros of / considered as a function of any of the Z\. 

The defining characteristic of a Laughlin wave function at filling fraction 1/q (where q is 
an odd integer and N s = qN e ) is that it vanishes as — Zj) q as any Z{ approaches any other 
Zj. When combined with the condition discussed in the last paragraph, namely that the 
sum of all zeros for each z^ should be constant, the only possible 
part of the wave function for the Laughlin state on the torus is 2], 



"orm of the holomorphic 





f(zi, z 2 ,...,z Ne ) = F cm ( ]T HiM^ - ^W- (8) 

i i<j 



Here F cm (Z) is a holomorphic function possessing q zeros, and is one of the four Jacobi 
theta functions: 

1/2 



1/2 

1/2 







1/2 



It), 



(9) 



where the theta function with characteristics is defined as 







r) = exp [i7rr(ri + a) 2 + 2ni(n + a){z + 6)]. 



(10) 



We will usually write the theta function with characteristics as $[a](£, r) where a is the 
vector (a, b) T . 

Each of the theta functions possesses a single zero in the period parallelogram, but note 
that for a, 6 G Z/2 

%](-z,r) = (-l) 4a6 %](^r), (11) 

so only $i obeys 2) = — and has its zero at z — 0. All theta functions are 
quasi-periodic, and, in particular, 

#i(z + l) = ~Mz), 

^(z + t) = -exp[-i7r(22 + r)]i?i(2). (12) 

Eqs. d3J), (|HD, and (|T2l) . together with iV s = giV, require the following quasi-periodicity 
conditions for F cm (Z): [2| 

F cm (Z + l) = (-lY N ^F cm (Z), 
F cm (Z + r) = (-lY N ^exp[-t7rq(2z + r)}F cm (Z). 

A convenient basis for such F cm (Z) is provided by the set of q functions 

m/q + (N s - q)/2q 
~{Ns ~ q)/2 



(qZ\qr), 



(13) 



(14) 
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where m is an integer defined mod q. The resulting set of q linearly-independent wave 
functions differ only by a rigid translation along the y direction, and in the limit N e — > 
oo each electron need move only infinitesimal amount to cause this shift. The distinct 
wave functions are locally indistinguishable, and local perturbations change the energy of 
these states by same amount. The g-fold degeneracy is therefore unaffected by such local 
perturbations. 

If now Nh quasiholes are inserted at Wi, then N s = qN e +Nh. Since / should now vanish to 
the first power as^j^ Wj, the only way to have the sum of the zeros of each Zi independent 



of {wj} is to set [2, lll| 



z 2 , . . . , z Ne ) = f£> (j2 * + ^p-j IlM* - *>i) II^(^ - z ^ q - ( 15 ) 

It is useful to include a purely {w^j-dependent part to the wave function normalization so 
that 

»<->(*, ■ • • , ■>.; R,, . . . , R K J = exp (- £ £j exp (- £ Jj) f£> + ^ 



Y[lMm - ™ 3 )] 1/s Y[M:, - wi) rj^ife - (16) 



X 

i<j i,j i<j 



where Wi = + ir]i)/L x and R = (£, rj). It was argued by Einarsson 14( that taking one 
quasihole around another on torus should result in the same phase factor on a torus as on 
plane and this factor is accounted for by rL<j[' i M u '« — w j)] l ^ q - The exponent of the new 
gaussian factor can be explained from the fact that the charge of a quasihole is e/q. These 
added factors in Eq. (fl6j) should therefore have converted all Berry phases into explicit 
monodromies. 

One can easily check that for {^}, boundary conditions of Eq. (II 6 1) is exactly that of 
Eq. ([I]), provided that N s of Eq. (fl3j) is equal to qN e + Nh- Under a straight-line analytic 
continuation z — ► z ± 1 we have 

[$i{z)} 1/q -> exp(±nr/q)[Mz)] 1/q (17) 

and as z — > z ± r we have 



[#1(2)] — > exp(=|=i7r/g) exp 



ITi 



— (±2z + t) 
1 
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Also note for the center-of-mass wave function: 



F<£\Z±l/q) = (-l)"« e ±*<("h-ff)/9 



exp 



m 

±2ni — 
Q 



F^(Z±r/q) = (-l) N *e ±niiNh - q)/q exp[-m(±2z + r/q^F^^iZ). 



(19) 



Taking a quasihole around one of the generators of the torus therefore effects the following 
transformations: under Rj — > Rj ± L x it we have 



j • • • ) *-iV e > 



r N e ; Ri , 



Rj 



exp 



m 2N h -q-l 

±2m I 1 

q 2q 



x 



(20) 



and under Rj — > Rj ± L y y we have 



- r at, ; Ri 



) • • • ) i iV e i '■W) • • • ) J - ,J 2; 



Rj 



_ ±in/q 



exp (=FiLy£i/qf 



x *( n * ±1 )(r lj ...,r iVe ;Ri,...,Ri,...). (21) 

Eqs. f TSU]) . (|21|) reproduce the operator transformation of Wen and Niu [_3] up to a gauge 
transformation: motion of the quasihole about one torus generator reproduces the state up 
io phase, and motion about the other generator rolls the ground state over into another one 



3JJ. These operations and their effects are analogous to the operations appearing in the 



Verlinde algebra of conformal field theory. 



III. MOORE-READ STATE 

The holomorphic part of the v — 1/2 Moore- Read state wave function on the plane is 

/mr(*1, ■ ■ ■ , ^2n) = Pf IB* ~ Z ^> ( 22 ) 

where the Pfaffian of a 2n-by-2n antisymmetric matrix is defined as 

1 ™ 

pf ( A ) = £ sgn(p) n^^po- ( 23 ) 

' P6S 2 „ »=i 

Here P runs over all permutation of 2n objects. This wave function is an antisymmetric 
polynomial in the Z{, the poles in the Pfaffian part having canceled against some of the 
zeros in the Laughlin-Jastrow factor. The suppression of these zeros by the Pfaffian factor 
increases the amplitude for the one electron to approach another. The Pfaffian can therefore 
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be considered to indicate a pairing between electrons [6j, and this pairing is closely related 
to the BCS pairing in p + ip superconductors [8|. 

In constructing the Moore- Read wave functions on a torus we must satisfy the boundary 
conditions Eq. ([3]). These boundary conditions come from the gauge choice and the fact 
that all electrons are in the lowest Landau level; they do not depend on what correlated 
many-body state the electrons are in. 

Now the Pfaffian part is equal to the correlator of 2n chiral Majorana fermion fields of 
the critical Ising model 

ty( Zl )...^))=Pf(—L-Y (24) 

V z * z j J 

There is more than one way to put such a correlator on a torus. This is because we are free 
to give the Ising fields periodic or antiperiodic boundary conditions around each generator. 
We will see that there is an intricate interplay between the boundary conditions ([3]) required 
of the physical electrons and the boundary condition choices we make for the Ising ^'s. 

A. Even-spin structure sector 

The Laughlin-Jastrow part of the wave function on the torus can be dealt with by the 
same (zi — Zj) — ■> $1(2^ — Zj) substitution as before. The l/fe — Zj) factors in the Pfaffian 
require a more careful treatment. Because the Pfaffian contains sums of products any z- 
dependent factors arising from the z^z + 1 or z^z + t quasi-periodicity properties will 
not appear as an overall common factor. Consequently we should arrange that there are no 
^-dependent factors at all, and simply setting l/(zi — Zj) — ► — z j) wn l n °t do. A 

possible choice is to set [6j 

I y[a](zi-Zj) ^ 



Zi Zj $l\Zi Zj j 

for some choice of a = (1/2, 0) T , (0, 0) T , (0, 1/2) T . The functions i?[a](z) / $i(z) now obey 
the following boundary conditions: 

#[a)(z + l) _ _ e2m J[a}(z) 



M z + t) M z 
These functions are proportional to the two-point functions 



WzWz')) (27) 
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for chiral Majorana fermions with different spin structures — i.e. with different periodic 
or antiperiodic boundary conditions round the generators of the torus. At least one of 
these boundary conditions must be antiperiodic. If this condition is not met, the fermion 
propagator has a zero mode, and the usual two-point function does not exist. We are 
therefore, at the moment considering only even spin structures, i.e. spin structures a = 
(a, b) T having Aab an even integer. It should be noted once again that these periodic and 
antiperiodic boundary conditions are merely properties of ingredients in the Pfaffian. We 
are not changing the boundary conditions of the many-body wave function. Instead the 
sign factors are accommodated by the center-of-mass wave function. This means that the 
center-of-mass wave functions obtained by setting q = 2 in Eq. (ITlj) . 



m/2 + (N s - 2)/4 
-(N s - 2)/2 



(2Z I 2t) 



(28) 



will no longer suffice. Instead there are distinct center-of-mass wave functions for each 
a = (a, b) T : 



m/2 + (N 8 - 2)/4 + (1 - 2a) /A 
-(N s - 2)/2 - (1 - 2b)/2 



(2Z\2t). 



(29) 



*<«."0(r lt . . . , r Ne ) = exp (- £ ^)F^(£ ^ Pf ( V^- t) ~ ^ 

V V l\ Z t Z j) J i^A 



Here m is an integer defined mod 2. (It will be shown later that Eq. (128]) is the center-of-mass 
wave function in the odd spin structure.) The complete wave function 

'•&[a]{zi- Zj)" 

i<j 

(30) 

with a = (1/2, 0) T , (0, 0) T , or (0, 1/2) T , now satisfies the boundary conditions Eq. (j3J). Note 
that the additional a-dependent terms that Eq. (|29|) posesses in comparison with Eq. (|28|) 
ensures that that the boundary conditions of the total wave functions are same for different 
a's. The three spin-structure choices coupled with the two possible values of m in Eq. (1301 
give us the sixfold degeneracy of the even-spin structure Moore-Read state on the torus 



IB 



io|. 



Now we consider inserting quasiholes. A charge e/2 quasihole excitation in one of these 
degenerate ground state is little different from a quasihole in a Laughlin state. Such a 
quasihole has one quantum of flux and consequently N s = 2N e + 1. When an e/2 quasihole 
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is inserted at w, the holomorphic part of the wave function becomes 



./ 



(a,m) . 



Zl, 



Z N e 'l 



i w J -'cm 



i i<j 

(31) 



When this quasihole is carried around the generators, the wave function Eq. (1311) behaves 
almost the same as the torus Laughlin wave function. Under w — > w ± 1 we have 



/ 



(a,m). 



Zl, 



z Ne \w) 



T^(-l) m f {a ' m) (z l ,...,z Ne ;w), a 







(32) 



and under w-t«)±rwe have 



/^)(^...,^ e ;«;) 



exp [-Hf^(±2 W + t)} f^ m+1 \z h ...,z Ne ;w), 



zpiexp [-Hf^(±2w + r)] /^ m+1 )(z 1 ,...,z^;u;), 



a 



a 




(33) 



This is to be expected because creating the e/2 quasihole breaks no pairs and so has no effect 
on the BCS pairing characterizing the Moore- Read state. (These transformations have been 
discussed by Oshikawa et al. in the operator language 10|. There, the monodromy was 
associated with the adiabatic insertion of a unit flux quantum into the "holes" of the torus.) 

The e/2 quasihole is not, however, the elementary excitation for the Moore- Read state. 
By allowing pair-breaking, a charge e/2 quasihole with one quantum flux can fractionalize 
into two charge-e/4 quasiholes, each possessing a half quantum of flux [4|, |8j. To construct 
a wave function with two charge e/4 quasiholes, the Pfaflian part needs to be modified, as 
this is the part that describes the BCS pairing. 

There are three basic conditions to be considered when writing down such a two-quasihole 
wave function. The first is that when these two quasiholes are brought together we should 
recover Eq. (I3"T|) . The second condition is that if two paired electrons wind around a single 
charge e/4 quasihole, it should result in accumulation of phase 2tt. The third condition is 
that, since the argument of the center-of-mass wave function becomes Zi + (w\ + u>2)/4, 
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the boundary condition of the Pfaffian part for the z% — > z% + r translation needs to change 
accordingly. Greiter et al. obtained wave functions satisfying these three conditions [6]. The 
holomorphic part of their wave function is 

f^ m \z u ...,z Ne] w u w 2 ) = Ft m) [J2 z * + Pf {M$) HMzi ~ z 3 )]\ (34) 

\ i / i<j 



with 

, , a _ $la](zi ~ Zj + wi 2 /2)#i(zj - w^Mzj - w 2 ) + (i <->■ j) , . 

M ^ ~ 2A(z i ~z j ) ' (35) 

and Wu = W\ — w 2 - (Note that (i <-> j) refers to a term that differs from the previous term 
only by the exchange of the index i and j; hence, for Eq. ( l35l) . (i <-> j) = i?[a](zj — Zi + 
Wx2/2)-&x{zj —Wx)'&i{zi — W2).) As far as translation of electrons around the torus generators 
are concerned Eq. (1341) has same transformation as Eq. ( 13TI) . or as the holomorphic part of 

Eq. dHDD- 

The above considerations are not, however, sufficient for obtaining the complete two- 
quasihole wave functions. Eq. fl34l) and Eq. fl35l) satisfy constraints involving only the electron 
co-ordinates 2^. There are also important factors that depend only on W\ and w 2 - To obtain 
the complete dependence on the quasihole coordinates as well, we need to calculate the 
holomorphic conformal blocks of the correlators for critical Ising model. 

The starting point is to observe from Moore and Read's original derivation how the 
Pfaffian part of the Moore-Read state wave function with two quasiholes is obtained when 



place on the plane {4]: 



Pf J - = {ip(zi) ■ ■ ■ ip{z Ne )a{wi)a{w 2 )) 



Z{ Zj 



X w 



H 8 l[[(z i -w 1 )(z i -w 2 )} 1 / 2 , (36) 



where ip is the Majorana fermion and a the spin field. One can immediately see that the 
Pfaffian part of the wave function on a torus would be equal to 

<V(*i) • ■ • ip(z Ne )a(^iH^2))a[A(wi2)] 1/8 HMzi ~ w 1 )A(z i - w 2 )} 1 ' 2 . (37) 

i 

Here a = (a, b) T denotes the boundary conditions, Eq. ([26]) . of the ip field on the torus, just 
as in the case of Eq. (J27j) . 
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The following Ising model correlators have been obtained in the even spin structures on 
the torus by Di Francesco et al. [ij], lfj 

(if;(z i )if;(z j )a(w 1 )a(w 2 )) a ^(0) 



(a(w 1 )a(w 2 )) a 2$ 1 (z i -z j ) 

/ \ 1/2 

•&[a](zi - Zj + W12/2) / tfi{Zi - w^d^Zj -w 2 y 



x 



^h|(>r 12 /2) \0v(:,,- u^'hlZj try) ' + " " 7 1 

(38) 



= to) ' 

(To be precise, these are the chiral holomorphic parts extracted from the non-chiral Ising 
field correlators obtained by Di Francesco et al.) From the antisymmetry of the ip field 
under exchange, and from the condtitions 



lim (zi - z 2 ) (if>(zi)if>(z 2 )if>(z 3 ) ■■ ■i){z Ne )a(w 1 )a(w 2 )) a = (ip(z 3 ) ■ ■ ■ip(z Ne )a(w 1 )a(w 2 )) a , 
lim ( Zi - ^(V'^W^^K^W). = 1; (39) 



Zi > Zj 



we can obtain the iV e -pomt if) field correlator (N e even) as 

mzi) ■ ■ ■^{z Ne )a[wi)a{w 2 )) a = {a{wi)a{w 2 ))aPf I (o-( Wl )a(w 2 )) / ' 

Eqs. fl37|) . fl38|) . and fj40|) tell us that the wave functions of Eq. (1M1) need to be modified in 
the following way if they are to have correct dependence on quasihole coordinates: 

f^ m \z 1 ,...,z Ne ;w 1 ,w 2 ) = F&H {^Zi + ^P^J [%](«W2)] 1/2 Pf (m?) 

i<j 



with 

jTra = 0[oi\{zi ~ zj + w 12 /2)$ 1 {z i - wQgife ~ w ^> + j) U2) 

2$ 1 {z i -z j )$[a]{w 12 /2) ■ 1 ' 

Note that while Eq. (|38|) is not single-valued in coordinates of the if) fields — taking a if) 

around a a results in a minus sign — the wavefunctiuon Eq. (|4ip is single-valued in the 

electron coordinate. This is because the [^(it^)] 1 ^ 8 ILI^M 2 ^ ~ w i)$i{ z i ~ w 2 )] l l 2 factor 

(in addition to the Laughlin-Jastrow factor) makes everything analytic in Eq. (EH), save 

for [ , i9[a](w;i2/2)] 1 / 2 . Thus the wave function is single- valued in electron coordinate, and 

has no singularity. Note that constants i?[a](0) and ^(0) have been ignored in Eq. ( )4Ti) . 



12 



(This means in the limit W\ — * w 2 Eq. (j3Tj) will differ from Eq. (13T|) by some multiplicative 
constant.) 

We now ask what happens when one of the quasiholes is translated around the generators. 
In contrast to Eq. ([3"3l where only m changed, this operation will result in a change in a. To 
see what happen s to the wave functions, we require the following standard theta function 



identities 
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1^: 



d 2 {z±l/2) = 

# 2 (z±t/2) = exp[-m(±z + T/4:)]$3(z), 

3 (*±l/2) = Mz), 

$ 3 (z±t/2) = exp[-i7r(±z + T/4)]$ 2 (z), 

ti 4 (z±l/2) = 3 (z), 

■d±{z±r/2) = ±iexp[-i7r(±z + r/4)]t? a (z), 

together with the center-of-mass wave function formulas, 



cm 


=2,m) 


± 1/4) 


_ e -i(l/2 T l/2)7r(m+7V s /2-l) p( m ) 




77.(0= 
cm 


=2,m) 


(z±t/4) 


= e ±™(N s -m exp[-m(±z + r/8)]F c ^= 


=3,771+1/2+1/2) f z 


77.(0= 
cm 


=3,m) 


(z ± 1/4) 


_ e -i(l/2 T l/2)7r(m+Af s /2-l/2) ^(a=4,m)/^ 




77.(0= 
cm 


=3,m) 


(z ± r/4) 


= e ±*^.-i)/4 exp[-z7r(±z + r/8)]F c ( m = 


=2,m+l/2±l/2) f z 


f( a = 
cm 


=4,m) 


(z ± 1/4) 


i(l/2±l/2)7r(m+Ar s /2-l/2) 77.(0=3,771) / _\ 




cm 


=4,m) 


(z±t/4) 


= e^^" 2 )/ 4 exp[-i7r(±z + t/8)]P£ 


•f-1/2+1/2)/ \ 



(43) 



(44) 

In these center-of-mass wave function formulas, the notation of Eq. (Q is used; that is, a = 2 
stands for a = (1/2, 0) T , a = 3 for a = (0, 0) T , and a = 4 for a = (0, 1/2) T . 

Note that the transformations of Eq. (1431) are covariant with those of Eq. (1441 . 
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Consider the results not involving $1 or F^': 

Wl ^ Wl ±r: f( a = 2 ' m \z 1 ,...,z Nts ;w 1 ,w 2 ) -> exp[-i7rN s (± Wl /2 + r/4)] 

x /( a = 3 ' m+1 /^ 1 / 2 )(^ 1 ,...,z^; Wl ,^ 2 ), 

^i-^±l:/ {a=3 ' m) (^,...,^ e ;^i,^) - (_l)^/2 e ^(i/2 T i/2) 

x /( a=4 ' m )( 2l ,...,^ e ; Wl , W2 ), 

w^^ir:/^^,...,^;^^) -> exp[-z7riV s (± Wl /2 + r/4)] 

x ^= 2 ' m+1 / 2±1 / 2 )(, 1 ,...,^; Wl , W2 ), 

mi -> w x ± 1 : / (a=4 ' m) (zi, . . . , zjv e ; Wlj m 2 ) -> (_i)^/2 e ^ m (i/2±i/2) 

x /( a = 3 ' m )(^ 1 ,...,z iVe ; Wl , W2 ). (45) 

In Moore and Read's original wave function formulation |4J the Gaussian factor comes en- 
tirely determined from the charge sector. Since the charge of a quasihole considered here is 
e/4, an analogy with Eq. (|T6|) shows that the total wave function should be 

¥ a > m \r u ...,r Ne ; R 1; R 2 ) = exp (- £ y 2 /2Z 2 ) exp [-( r/l 2 +r / 2 )/8/ 2 ] /<«•"»>(*!, ■ • • , z N .\w ly w 2 ). 

(46) 

Eq. (I45p and Eq. (j4"6l) gives the following transformation rules for the total wave functions: 

R x ^ R x ± L y y : ^= 2 < m )( ri , . . . , r Ne ; R 1; R 2 ) - exp ( T iL y ^/4l 2 ) 

x ^ a=3 ' m+1 /^ 1 / 2 )(r 1 ,...,r iVe ;R 1 ,R 2 ), 
Ri -> Ri ± L..X : *( a = 3 ' m )(r 1 , . . . , ix; Ri, R 2 ) -> (_i)^/2^ m (i/2 T i/ 2 ) 

R x - R x ± L y y : ¥ a=3 > m \ ri , . . . , r Ne ; R 1; R 2 ) - exp (^L^/4/ 2 ) 

x ^= 2 ^ +1 / 2±1 / 2 )(r 1 ,...,r 7Ve ;R 1 ,R 2 ), 

Ri -> Ri ± L,.x : $( a=4 ' m )(n, . . . , ix; Ri, R 2 ) = (_i)^/2 e ^m(i/2±i/2) 

x ¥ a=3 ' m \r 1 ,...,r Ne ;R 1 ,R 2 ). (47) 

We see that, unlike for the Laughlin states, these states are not eigenstates with respect 
to Wj Wj ± 1 transformations. However, the feature that this transformation does not 
change m persists. For the translations Wj Wj ± r, there is a phase factor due to the 
gauge transformation A — > A — BL y it. Comparing the gauge transformation phase factors 
of Eq. (1211) to that of Eq. (j4"T|) confirms that the quasiholes in Eq. (J4"T]) have charge e/4. 
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Eq. (l4"5j) left the following two cases unmentioned: 

W 1 ^W 1 ±1: f {a=2 > m \ Zl , . . . , Z Ne ; Wl , W 2 ) -> (_l)^/2 e --(l/4 T l/4) ei7 r(m-l)(l/2 T l/2) 

x / (m) (2; 1 ,...,^ e ;M;i,w 2 ), 
wx^^ir:/^ 4 ^^,...,^;^,^) - e^/ 4 wp[-mN s {±w x /<2 + r/4)] 

x /( m+1 / 2±1 / 2 )(^ 1 ,...,^;«; 1 , W2 ), (48) 



where 

f^(z u ..., ZN ,, Wl , W2 ) = s [^(W2)] i/2 Pf (-Mr (1/2,1/2) ) 

x n^(^-^)] 2 - (49) 

The tilde mark is placed above / in Eq. P9]) because, in some sense, it does not qualify as a 
Moore- Read state wave function. When the two quasiholes are merged in Eq. ( H9|) . the wave 
function simply vanishes; one does not get one of the wave functions of Eq. (13T1) . This does 
not mean that Eq. (1491) necessarily gives higher energy than Eq. fj4~T]) . However, note that 
Eq. fl31|) differs from the ground states only by one quantum flux; Eq. fT49|) does not have a 
corresponding ground state in this sense, hence its disqualification. One can further regard 
having two e/4 quasiholes as differing from having a quasihole-quasiparticle pair merely by 
one quantum flux. This leads to the conclusion that if one creates a quasihole-quasiparticle 
pair out of one of a=2 (or a=4) ground states, translate the quasihole around the generator 
in x (or y) direction, and annihilate the pair, one does not return to a ground state; this 
means that the monodromy process can actually excite the system. 

To find out what kind of excitation do we have here, we need to examine the vanishing of 
Eq. ( 1491) as two quasiholes are brought together. An examination into Eq. (1421) shows that 
the vanishing is not due to the Pfaffian part of Eq. ( 1491) . Indeed 

lim M a=< > 1/2 ' 1/2) - -wijtiijzj -wi) \$'i{zi - Zj) ( tf'^Zj - wt) ^[(zj - wi 



- 2 —i 13 #i(0) W*-u>i) M*i-vi) 

(50) 

The vanishing of Eq. (141)1) is therefore solely due to the [$i(wi 2 /2)] 1//2 . This term originates 
from the correlator of two quasiholes in the a = (1/2, 1/2) T sector. This is the sector with 
the periodic boundary conditions around both generators for ip fields: 

(«r(wi)<r(«> 2 )) a ^ (1/2 ,i /2) [i?i(«Ji 2 )] 1 / 8 oc [0i(W2)] 1/2 - (51) 
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The spin structure is therefore correlated with the internal state (fusion channel) of the two 
quasiholes. Since for w — > 



(<7H<7(0)) a=(1/2 , 1/2) ~ w 3/8 , (52) 

whereas for a = (1/2, 0) T , (0, 0) T , and (0, 1/2) T 

(a(w)a(0)) a ~ w- 1 / 8 , (53) 

the chiral Ising model operator product expansion 

a(w)a(0) ~ — ^ + const.w 3/ V(0) (54) 

tells us that for a = (1/2, 1/2) T two cr's fuse into r/>, whereas for a = (1/2, 0) T , (0, 0) T , 
and (0, 1/2) T they fuse into I. We see that our wave function approach makes manifest the 
observation of Oshikawa et al. that, after the translation of a quasihole around a generator, 
the system ma y r efuse to go back into a ground state because of a change in the a x o 



fusion channel 



10] . This change was explained by Oshikawa et al. by using the branch cut 



argument of Ivanov 17| and Stern et al. 181 ] 



From the preceding arguments, one can explain this change of the fusion channel in the 
terms of the physical fusion in the manner discussed by Stone and Chung for the two- 
dimensional p + ip superconductor [3]. In the ground states, two cr's fused into I and the 
total fermion number on torus was even, enabling all fermions to be paired up. However, 
once the translation of one quasihole rolls the system over to the a = (1/2, 1/2) T sector, the 
cr's fused into a ?/>, which, due to the conservation of total fermion number, is possible only 
if there is depairing of one of fermion pairs. Therefore, with this change of fusion channel, 
the system is left in an excited state. 

This superconductor analogy also indicates that the parity of electron number in the 
ground state of the a = (1/2, 1/2) T sector should be different from that of other sector. In 
this picture, after the roll-over to the a = (1/2, 1/2) T sector, there is a Bogoliubov quasi- 
particle excitation in the system. Now if superconductor with one Bogoliubov quasiparticle 
excitation has iV electrons, this is equivalent to having a superposition of one hole excitation 
on a ground state with N + 1 electrons and one particle excitation on a ground state with 
N — 1 electrons. Thus the change of fusion channel has to be accompanied by the change 
in the parity of the electron number in the ground state. 
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As previously commented, in discussing topological features, creating a quasiparticle- 
quasihole pair is equivalent to splitting a charge e/2 quasihole into two e/4 quasiholes; the 
only difference between two cases is that there is one more flux quantum for the latter case. 
In their paper, Oshikawa et al. defined T XtV to be a process in which creates a quasiparticle- 
quasihole pair is created and then the quasiparticle is dragged around the generator of the 
torus in x (or y) direction to wrap around the system before it is pair- annihilated with the 
quasihole 10j. So one can make following correspondences between the analytic continuation 
of theta functions considered in this paper and the processes defined by Oshikawa et al.: 

wi -> wi + 1 T x , 

wi — > w\ + r <^ T y . (55) 

It will be shown in the appendix that the ground states labeled by a and m in this paper are 
eigenstates of T~ 2 and Ty. Labeling eigenvalues for these operators divided by the gauge 
transformation factor as f y and f' x respectively, we obtain the following correspondences 
between the wave functions of this paper and the state vectors of Oshikawa et al.: 

* (a=2 ' m=0) <- \fy =i,f x = 1) tf(«=2,-=D <_> \f y = _i, // = 1), 

^(a=3, m =0) ^ ^ = h f 9 = 1) y(a=3,m=l) „ ^ = _^ = ^ 

(-1)^/^=^=0) „ lfy = 1, // = -1) ( _ l} iVe/2^(a=4, m =l) ^ ^ = _ Xj = 

(56) 







1/, = 1, /* 


= 1) 


%c\fy 


hf' x 


= 1) 


= \fy 


%j\fy ~ 1> 


£ = 1> 


= |/» = -i, / 


:. = i) 


T X Ty\fy = 


hfx 


= 1) 


= o, 




/:=i> 


= \f v = -l,J 


; = i) 


%? y 2 \fy = 




= 1) 


= 1/, 


Tylfy = 1, 


= 


= \fy = h fx 


= 1) 


Wy\fy = 




= 1) 


= 0. 



Adapting the continuation formulas of Eq. (|48|) to the correspondence made in Eq. (job)) , we 
find the following topological actions 

1, /£ = -!>, 
-1, = -1>, 

(57) 

There are exactly the transformation formulae of Oshikawa et al. [10] Note, however, that for 
our formula, there is a provisio that at after all the quasihole translation has been carried 
out, the two quasiholes are to be merged . As in Eq. (j2T)) . we have dropped the gauge 
transformation phase factor. 
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Oshikawa et al. considered bases diagonalized with respect to either T x or T y [lCj . These 
two bases are related by the modular S'-matrix. In either basis, it can be shown that the 
transformation formula? of this subsection follow the Verlinde formula on diagonalization of 



fusion numbers by the S'-matrix 
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B. Odd-spin structure sector 

We now need to consider the odd spin structure in detail. On a torus, this is the case where 
the chiral Majorana fermion has periodic boundary conditions around the both generators. 
The action of a chiral Majorana fermion field ip is 



S = — d l zi)di\). (58) 
2n J 

In performing the path integral, it is necessary to expand ip in terms of the normalized 
c-number eigenmodes ip nm (z, z) of 8,. These must be doubly periodic, and so are 

1 71 

i> nm {z, z) = exp [n(rz - tz) + m(z - z)\. (59) 

V Im r Im r 

The resulting mode expansion is 

if)(z, z) = ^ a nm ip mn (z, z) = a 0Q + ^ \a nm ^nm{z, z) + a (z,z)}, (60) 

n,m£Z n,m>0 

where a n 's are Grassman variables with a_ n _ m = a* nm and Yl'nm ^ s a summation over non- 
negative integers that excludes the n = m = term. It should be noted that since if) is 
holomorphic only as a result of the equation of motion, one cannot take it to be holomorphic 
in performing the path integral. 

When the mode expansion Eq. (|60|) is inserted into the action of Eq. (|58|) . we obtain 

1 

S = — } (nr- m)a- n m a nm . (61) 

n,m>0 

The Grassmann variable aoo does not appear in this action. However, aoo does appear in 
the integration measure 

d[V>] = da nm . (62) 

n,m£Z | 1 | 1 | 1 

This leads to the odd spin-sector partition function vanishing 15j, [l6|, l2l| : 



/ da 00 I | [ da nm da_ n _ m J exp(-S') = 0. (63) 

\n,m>0 J 
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As in the Yl'nm °f ^q. dHO]) , Y[' nm i s a product of terms with non-negative integers n, m that 
excludes the n = m = term. Now J a da = 1 and J da = for any Grassmann variable a, 
and so the integration J F(a) Ylnm^ a nm yields zero unless each a nm appears exactly once 
in the integrand F(a), and there is no aoo in the partition function integrand. For the same 
reason, any correlator of an even number of ^'s vanishes. 

Correlators with odd number of ?/>'s can be non-zero. Consider the simplest case: 



{ip{z)) oc / d[ip]ijj(z) exp(— S) oc 



/ a 00 da 00 I J~J da nm da_ n _ m J exp(-S') 

\n,m>Q ) 
—n,—m 

exp(-S) ^ 0. 

\ n,m / 



(64) 



Thus (ip(z)) is a nonzero constant [15I, Q]- The next simplest case is the three-point corre- 
lators, where three ip fields take turns in occupying the zero mode: 

(ip(z 1 )ip(z2)ip(z 3 )) oc / d[if>]iJ)(z 1 )if>(z2)il>(z 3 )exp(-S) 



oc / ( Yl da nm da_ n _ m \ ^ a^tp^^)^ a M ip k i(z 2 ) exp(-S) 

J \n,m>0 / i,j>0 k,l>0 

+ / ( II da nmda_ n ^ m ) ^ a ij ij ij (z 2 )^2 a k iij k i(z 3 )exp(-S) 

J \n,m>0 / i,j>0 k,l>0 

+ J ( II da nmda_ n _ m ) ^2 a ij^ij( z 3)^2 a kii'ki(z 1 )exp(-S) 



\n,m>0 / i,j>0 k,l>0 

g ++ (r l - r 2 ) + g++(r 2 - r 3 ) + g++(r 3 - n), 



where 



BE 



9++{ r ) 



A(z) L y ■ 



is a modified Green function satisfying 



71 

Imr 



(65) 



(66) 



(67) 



This Green function contains y = (z + z)/2, and is not holomorphic. The f's, however, 
cancel in sum of three Green functions appearing in Eq. ( |65l) . The 3-point correlator is, 
therefore, holomorphic. 

In the odd spin structure, and for N e odd, the general iV e -point chiral Majorana fermion 
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correlator is 

(V>(*i)...V(*jO) 00 ■■■>*jO 



1 ^ (iV e -l)/2 

2(j v e -i)/2 [(iVe _ 1)/2]! E s g R ( P ) II ^++(rp(2 4 -i)-r P(2i) ). 



(68) 

Here P runs over all permutations of N e objects. Again, at first sight,this equation looks 
neither chiral nor holomorphic. However, Read and his coworkers 3, |8j asserted that the 
cancelation in Eq. fl65l) generalizes to 

...,,) - „^ _ mi x «■<*> Y i£H : £Si • (69) 

This rather non-obvious cancelation of the z's will be proved in the appendix. 

Eq. (j68p was first constructed by Read and Green in their work on spinless p + ip super- 
conductor 8]. They have pointed out that the existence of the k = mode, which is an 
almost exact analog of the zero mode discussed here, requires the number of electrons to be 
odd, as long as it is energetically favorable to have this k = mode occupied. 

We emphasize again that these periodic and anti-periodic "boundary conditions" for the 
^(z)'s are not the boundary conditions of the physical electrons. The physical boundary 
conditions are those of Eq. (JTj) and The boundary conditions for the ^(^)'s affect 

only the manner of BCS pairing. The best way of regarding Eq. ( 168]) is that this is the 
correct replacement for the Pfaffian when the number of electrons is odd. The center-of- 
mass wave functions ensure that all Moore-Read states on a torus have the same physical 
boundary conditions under translation of electrons around generators, regardless of whether 
the number of electrons is even or odd. This means that F c ( m } of Eq. ([28]) is the center-of- 
mass wave function when the number of electrons on the torus is odd. The holomorphic 
part of the ground state wave function in this case is therefore 

f!^(zi,...,z Ne ) = 7-r f >> ) -v. ) I |KU, - :,)] 2 . (70) 



( e * ) ■ • • > n^^ - z ^ 

\ i / i<j 



Obtaining the wave functions with one charge e/2 quasihole with one quantum flux is now 
straightforward. It is in essentially same form as Eq. ( |3T]) : 



\ i / i i<j 



(71) 
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How one do we fractionalize the quasihole of Eq. (1711) to obtain two charge e/4 quasiholes? 
The first formula of Eq. (139|) . being independent of the parity of number of t/>'s, still holds. 



Also note if we define M° dd 



lim 



-M. 



a=(l/2,l/2) 



Kl 



then 



[zi-zM(0)MS 



odd 



1. 



(72) 



[$i{zi - w^d^Zi - w 2 )^i{zj - w^tftizj - u; 2 )] 1/2 
Therefore, in the odd spin structure for N e odd, the chiral Ising field correlator of Eq. (|40| 
becomes 

(ip{zi) ■ ■ ■ ip(z Ne )a(wi)a(w 2 )) a =(i/2,i/2) = - - w 2 )]~ i/2 

8=1 

P&S Ne 

x (4>(z P ( Ne ))cr(wi)cr(w2))a=(i/2,i/2) 

(N e -l)/2 



X 



n 



iW P(2i-l),P(2j) 



(73) 



Note that one now needs to calculate (^o"c r ) a =(i/2,i/2) rather than (<To) a =a/2,i/2)- This leads 
to the conclusion that with two e/4 quasiholes at Wi and w 2 , Eq. (ITT!) should be modified 
into 



&++(zi, ■ ■ ■ ,z Ne ;w 1 ,w 2 ) 



const. sgn(P)[?9i(wi — w 2 )] 1 ^ 8 

Pes Ne 

x [^(^ e) -^ 1 )] 1 / 2 [^(^ (iVe )- W2 )] 1 / 2 

(AT e -l)/2 

X (i>(zp(N e ))v(wi)(T(w2)) a=W2 ,l/2) J| tet2i-l),P(2i) 



i=l 



We compute the necessary Ising correlator in the appendix. We find that 

1 1 1/8 

(tfj(z)a(w l )a(w2)} a =(i/2,i/2) oc — 

IV 1 {W12)_ 

n > / / n In/ / n {$1 (z — Wi) $1 (z — Wo) 

$ x {w l2 2) + -^K 2 -i) ^ - -i) ^ 

2 \vi{z — wi) Vi(z — w 2 ) 



Consequently, ignoring a multiplicative constant, we have 

(iV e -l)/2 

^ ++ {z 1 ,...,z Ne ;w 1 ,w 2 ) = ^ sgn{P)[h{z P{Ne y,w 1 ,w 2 )} 1 / 2 JJ [M P d { i_ 1):P{2l) 

PGS Ne 1=1 



(74) 



1/2 



(75) 



(76) 
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where 



h(z;wi,w 2 ) = $' 1 {w 12 /2)<& 1 {z-w 1 )'& 1 {z-w 2 ) 

+ h x {w X2 /2) (d' 1 {z-w l )d l {z-w 2 )-d' 1 {z-w 2 )d 1 {z- „-,)). (77) 

The holomorphic part of the wave function with two charge e/4 quasiholes is 

<.(m)/ s 



\ i J i<j 

\ i / i<j 

J2 sgn(P)[h(z P{Ne y, Wl ,w 2 )}^ J] [M° P d $ k _ 1)>P{2k) ] . (7S 



(iV e -l)/2 



2 



X 

PeSjv, fc=i 



It is not obvious that this wave function Eq. (1751) is analytic in electron coordinates. We 
need to show that [h(z; Wi, w 2 )\ 1 ' 2 is an analytic function of z. The proof of the analyticity 
of [h(z; w\,w 2 )} l l 2 can be presented in two steps. First, note that from 

h(z + 1; w\,w 2 ) = h(z;w 1 ,w 2 ) 

h(z + t; wi, w 2 ) = exp[— 2iri(2z — wi — w 2 + r))h(z; w%, w 2 ) (79) 

the z of h(z; W\,w 2 ) should have 2 zeros in the first principal region. Second, note that both 
h(z;wi,w 2 ) and d z h(z;wx,w 2 ) vanish at z = (wi + w 2 )/2. This shows that there is actually 
a double zero at z — (w\ + w 2 )/2 and this is enough to ensure that [h(z; w\, w 2 )] l l 2 is indeed 
analytic in z. 

We should also verify that the odd-spin-sector wave function with two e/4 quasi- 
holes goes smoothly into the corresponding wave function with one e/2 quasihole (up to 
a overall multiplicative constant) as the e/4 quasiholes merge. In other words, in the 
limit that Wi,w 2 — > w, does ifr ++ (zx, . . . , Zn b ; w%, w 2 ) of Eq. (1761) become proportional to 
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g++(zi, ■ ■ -,z Ne ) Ui^i( z i ~ W V- From Eq. ([SO]), 

(JV e -l)/2 



hm^+^i,...,^^,^ = moWUMti-v) E s S n ( P ) II 

,W2 ^ W i Pes Ne i=i 1(1 ' 

^i(^P(2i-l) - Z P{2 i)) ( fi'tiZp^i-l) ~ W) ^[(zpQi) - W)' 



' #l(Zp(2i-l) - 2P(M)) V ^(^(W-l) ~W) l?l(^P(M) -W). 

(AT e -l)/2 Q/ , , 

= const. n^-« ) e ^) n t 1 ^-"'^ 

= COTJSt.'9 ++ (zi,...,Z Ne )Yl'&i(Zi-w). (80) 

i 

The cancelation of the w-dependent terms other than Y[ t "&x{zi — w) occurs for exactly same 
reason as the one in Eq. fl69|) . This cancelation will be proved in the appendix. 
From 

h(z;w 1 ±l,w 2 ) = ±$' 2 (w 12 /2)'d 1 (z — Wi^^z - w 2 ) 

± ^2(^12/2) K(z - Wl )A(z - w 2 ) - ^(z - w 2 )A{z - Wl)) , 
h(z;w 1 ±r,w 2 ) = ±ie' iTT/i e TiTWl2/2 e- inT e ±2niz [^' A (w 12 /2)^ l (z - w x )#i(z - io 2 ) 

+ - - » 2 ) - ^'1(2 - ^2)^1(2 - (81) 

it is clear that after after one quasihole is translated around a generator, the wave function 
vanishes when two quasiholes are brought together. 

As in the even spin structure case, this vanishing can be attributed to the change in the 
fusion channel of two cr's. This is so since 



lim (ip(z)a(wi)a(w 2 )) a=(1/2A/2) ~ w^ /8 , (82) 



Wi— VW2 



but after the translation W\ — > w\ ± 1 or w\ — > u>i ± r, Eq. (IHTl) tells us that the same 
correlator vanishes as w\!^ when W\ — > w 2 . From the chiral Ising operator product expansion, 
Eq. fl5^|) . one can see two cr's fuse to I in (ipcrcr) a =(i/2,i/2)-, but after the W\ — > W\ ± 1 or 
Wi — > iyi±r translation, they fuse to ■0. The argument made in the last subsection that this 
change in the fusion channel is accompanied by the change in the parity of the number of 
the electrons in the ground state remain valid. Since the number of electrons in the system 
does not change, a quasiparticle excitation must have been created. Note that the different 
monodromy results of the odd spin structure arises from the fact that here (ipaa) plays the 



23 



role of {era) in the even spin structures. We see that the different monodromy outcomes 
originate from the different parities of the electron number. 

From the monodromy outcomes we have obtained, the ground states of the odd spin 
structure are eigenstates of both T x and T y . However, in the even spin structures, no ground 
state is simultaneously an eigenstate of both T x and T y . The difference comes from the fact 
that, in the language of the critical Ising model, the total topological charge of the system 
is if) in the odd spin structure whereas it is I in the even spin structures. That means that 
the modular S"-matrix of the system for the odd spin structure differs from that of the even 
spin structures - for the former and S 1 for the latter 22J, |23|] . 

This discussion on total topological charge indicates that the odd spin structure sector 
exists only because electrons that makes up a quantum Hall system are fermions. On the 



other hand, a system consisting of bosons, as in the string- net of Levin and Wen 24| , one 
cannot have a ground state with total topological charge of if). In such system, one would 
still have 'forbidden transition' due to change in the fusion channel, but this change is not 
tied with the change in the parity of number of particles in the ground state in the manner 
discussed here. The change in the fusion channel can be attributed to the change in the 
parity of number of particles in the ground state only when the particles are fermions and 
the ground state has BCS pairing. 



IV. CONCLUSION AND DISCUSSION 



In this paper, we constructed the two-quasihole wave functions for the Moore- Read quan- 
tum Hall state with on a torus for both even and odd spin structures. Conformal field theory 
calculations enable us to obtain the complete dependence of these wave functions on quasi- 
hole coordinates. We showed that the number of electrons in a ground state must be odd in 
the odd spin structure, and even in the even spin structures. By noting that the boundary 
conditions of each electron remain that of Eq. (TjQ) for all topologically degenerate ground 
states, we obtained the explicit expressions for the center-of-mass wave functions in all cases. 
Analytic continuation allowed us to obtain the monodromy matrix describing the effect on 
the space degenerate ground states of taking quasiholes around the torus generators. The 
effects are in agreement with those obtained by Oshikawa et al. In this process, we demon- 
strated that otherwise anticipated transitions are forbidden because they involve a change 
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in the fusion channel of the quasiholes. This reflects a change in the parity of the number of 
electrons that can reside in the ground state. Since the number of electrons is conserved, the 
operations that might have resulted in these forbidden operations take us out of the space 
of degenerate ground states, and into the space of excited states. 



Several extensions of our ana^ 



ysis are possible. Theta functions can be generalized to 



higher genus Riemann surface |25j , and these functions will naturally be ingredients of wave 



functions in such topology. Such wave functions for Laughlin states had been studied [26]. 
With these wave functions are found for the Moore-Read state, it should be possible to 
extend our calculations to the higher genus Riemann surface considered by Oshikawa et 



al. 10]. It would also be valuable to extend this analysis to consider the next simplest 
non-Abelian quantum Hall states, the Read-Rezayi parafermion states. The number of 
degenerate ground states has already been found 27|. The challenge is to work out the 



wave function on the compact Riemann surface using the conformal field theory analysis, 
thus extending the recent wave function construction by E. Ardonne and K. Schoutens on 



the plane [28[] 
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APPENDIX A: TRANSLATION EIGENVALUES OF WAVE FUNCTIONS 

Following Eq. (}55l) . T~ 2 and Ty can be regarded as implementing Ri — > Ri — L x x 
twice and Ri — > Ri + L y y four times, respectively, on two-quasihole wave functions 
$' a,m '(ri, . . . , rjv e ; Ri, R2)- For a = 3, repeated application of Eq. (Ho^) leads to 

T y A :vl/ (a=2 ' m) (ri,...,r i v e ;Ri,R2) - exp(-zL^i// 2 )^ (a=2 ' m) (ri, • • • , r Ne ; R 1} R 2 ), 

r x - 2 :^= 3 .-)(r 1) ...,r JVe ;R 1 ,R 2 ) - e^ m ¥ a = 3 ^( ri , . . . ,r Ne ;R u R 2 ), 

T 4 : ^=3,m) (ri) ___ )rjVe . Ri)R2) ^ exp(-zL^ 1 // 2 )^( a=3 ' m )(r 1 ,...,r^;R 1 ,R 2 ), 

r; 2 :tM,™)( ri) ... )r%;Rl)R2 ) ^ e l7rm vl/(«=4,-)( ri) ... )r7Ve;Rl)R2 ). (Al) 

To apply these operations on a = 2 and a = 4 states, one also needs to consider the 
transformation of p m \z\, . . . , Zn b ', Wi, w 2 ) of Eq. fj4"§j) . even if this does not qualify as a 
Moore-Read state wave function. From the standard theta function identities 

0i(*-l/2) = -M*), 

•& x {z + t/2) = iexp[-m{z + r/A)]^ 4 (z), (A2) 
and the center-of-mass wave function transformation 
F {jn \z — 1/4) = F c ^= 2 ' m )(z), 

F^ m \z + r/A) = e in{Ns - 2)/4 ex P [-in(z + r/8)}F^ m \z), (A3) 

one obtains 

w x ^ Wl -l: f( m \ Zl , ...,z Ne ; w u w 2 ) - (-1) N ^ 2 f^ m \z u ...,z Ne ; w u w 2 ), 
w 1 ^w 1 + r:f^\z 1 ,...,z Ne ;w u w 2 ) -> e—/ 4 exp[-i7riV s ( Wl /2 + r/4)] 

x J^ 4 ^,...,^;^,^). (A4) 

Repeated application of Eqs. (j4"5|) and (1A4|) leads to 

^-2 . ^(a=2, m)(ri; ; rjVe . Ra) ^ _ ie i ff(TO -l)^(«=2,m) (rij . . . , rjv e 5 Rl, R2), 

T y 4 :^ a=4 - m )(r 1) ...,r JVe ;R 1 ,R 2 ) -> - exp(-i J L y a/^)* (o=4,m) (ri, • • • , rjv e ; Ri, R 2 ). 

(A5) 
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APPENDIX B: PROOF FOR CANCELATION IN SUMMATION 



Consider an function G(r, r') odd under exchange of r and r'. Suppose this function can 
be expressed as a sum of two function g x (r, r') and <? 2 ( r > r '); 9z{ v i v ') nas a further property 
that (yf 2 ( r i, r 2) +5'2( r 2, r 3) +fl , 2( r 3? r i) = 0. Certain functions in Sec. Ill, such as g++(r) in 
Eq. (1681) and M° dd in the limit w±, W2 — > w as in Eq. (|80|) (when regarded as a function of 
two electron coordinates), are of this form. In order that the assertions made in Eq. (168|) 
and Eq. (IHUj) hold, it is necessary to show that for N odd 

(iV-l)/2 (N-X)/2 

£ sgn(P) JJ C^p^-i),^)) = £ sgn(P) JJ &i(r P (2i_i),rp( 2 j)). (Bl) 
PeS N i=l Pe5]v i=l 

One can easily see that it should be so for N = 3 from G(rx, r 2 ) + G(r 2 , r 3 ) + G(r 3 , ri) = 
giiji, r 2 ) + gi(r 2 , r 3 )+gi(r 3 , ri), due to the property of g 2 mentioned above. Suppose it holds 
for N = 2k + 1 where /c is some positive integer. If it can be shown from this assumption 
that Eq. (1B1[) holds for N = 2k + 3, then we have a proof by induction. 

To apply induction, a new permutation P' G S 2 k+i needs be introduced: 
fc+i 

£ sgn(P) n G ( r ^-i)>rp( 2l )) = 2£(-l) m -"- 1 G(r m ,r n ) 

x £ sgn(P')n G ( r ^( 2 - 1 )' r ^'( 2 0)' ( B2 ) 

P'e5 2fc+1 i=i 

where m, n are integers between 1 and 2k + 3. This P' is a permutation of integers between 
1 and 2k + 3 except m and n. Since the assumption had been made that Eq. fIBll) holds for 
iV = 2A; + 1, Eq. (JB2]) means 

fc+i 

sgn(P)J]G(r P(2i _ 1) ,r P(2i) ) = 2 £ (-1)— ^(i^, r n ) 

PG5 2 fe+3 i=l m<ra 

fc 

x £ sgn(P') JJflfi(rp/( M _i),rp/(ai)). (B3) 
P'eS 2fc +i »=1 

Comparing Eq. flB2j) and Eq. (1B3|) shows that Eq. (1B3|) means we can replace all but one G 
into g\. However, this leads to the conclusion that the last remaining G can also be replaced 
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by gr- 

fc+1 fc+1 

^2 sgn(P) JJC(rp(2i-i),rp(2 i )) = sgn(P)G(r P(2fc+1) , r P(2fc+2) ) JJ #i(r P ( 2 i-i), r P(2 ^ 

P£S 2fe+3 i=l PeS 2k +3 »=2 



- sgn(P)[G(r P(2fc+1) , r P(2fc+2) ) + G(r P(2fc+2) , r P(2fc+3) ) 



3 

P£S 2k+3 



fc+1 



+ G(r P ( 2fc+3 ), r P(2fe+ i))] Yl 9i{rp(2i-i), r P ( 2 *)) 

= \ s g n ( P )M r P(2fc+i), r P(2fc+2)) + #i(r P(2fc+2) ,r P(2fc+3) ) 

P 652^+3 

fc+1 

+ 5 , l( r P(2fc+3), r P(2fc+l))] J|fl , i(rp(2i-i),rp(2i)) 

i=2 

fc+1 

= s S n ( P )Yl9i(r Pi 2i-i),r Pi 2i))- (B4) 

Pe5 2fc+3 i=i 

APPENDIX C: CRITICAL ISING CORRELATORS ON TORUS 

To calculate 

(i>(z)a(w 1 )a(w 2 )) a =(i/2,i/2), (CI) 

first note the connection between the correlator in the chiral Ising theory and the full Ising 
theory: 

\{ip(z)a(w 1 )a(w 2 ))\ 2 = (e(z,z)a(w 1 ,H) 1 )a(w2,w 2 )). (C2) 

So (ipacr) can be computed by taking the holomorphic part of (eaa). By the Ising model 
bosonization formula flol I30I]: 

1 ( -\ ( - w - \\i Q/<w 0(w 2 ,w 2 ), 
z)a{w 1 ,Wi)a{W2,w 2 )) = -2{d(p{z)0(p{z) cos cos ) 

= -^[(d4>(z)d(j)(z) expi^i^W!, w 1 )/2) exj>(-i(f>(w 2 ,w 2 )/2)) 
+ (d(j)(z)d(j)(z) exp(-z0(w 1; w x )/2) exp(i0(w 2 , w 2 )/2))], 

(C3) 

where z) is a free boson field. It should be emphasized that there is an intricacy hidden 
behind Eq. (IC3p . Critical Ising fields on torus can be bosonized only into a compactified free 
boson. A boson with compactification radius r = 1 should have boundary condition 

<j)( m > m '\z + l,z + l) = <j) m ' m '{z,z) + 27cm, <f) {m ' m/) {z + T,z + f) = f m '(z,z) + 2W, (C4) 
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where m, m' are integers. This winding comes entirely from the zero mode, so one can 
express the free boson field as a sum of this zero mode and "free part" - that is, nonzero 
modes: 

K ' ' ' (Co) 



^ m > m '\z, z) = —[m(rz - fz) + m'(z - z)\ + <p(z, z), 
where r = T\ + ir 2 with t\,t 2 real. (Note that in this appendix, the first principal region is a 



parallelogram and not necessarily a square and nonzero T\ is considered.) 
of the boson field leads to the following decomposition of the action: 



15 



lis decomposition 



2l| 



S[<j>] = (l/Srr) J (dfiW) 

= S[</> ] + S[$\ - (1/4tt) J 0A0 O = S[<f> ] + S[$\, (C6) 

where 0o refers to the zero mode terms of Eq. (j05j) . This result is due to the vanishing of 
the Laplacian of 0o- This also means that the boson partition function Z bos factorizes into 
the zero mode part Zq and the the "free part" Z: Z bos = ZqZ. Since in calculating Zq, all 



15 



16 



m, m! need to be summed over, one obtains 

ry \ ^ 7 (m,m') 

z q ~ 2 / ° 

m,m' 

where Z < ^ n,m ^ = exp[— (7r|mr — m'| 2 /2r 2 )]. 



2l| 



(C7) 



To compute a full boson correlator, one should be calculate the correlator first for the 
winding sector m,m' and then sum over all possible winding numbers. However, this does 
not suffice for the purpose here, which is to calculate (ipaa) Q =(i/2,i/2)- The question is how 
to extract out the portion of the correlator that corresponds to a = (1/2, 1/2) sector of the 
Ising model. 



The answer can be obtained from converting 



expressing theta functions in terms of Zq 



(m,m r ) 



of Eq. (lC7p . Following Di Francesco et ai, 16] 

2 





1 






2 


{z\r) 




1 






2 





e j7r[r(n+l/2) 2 -f(n+l/2) 2 ] e 27rj[(n+l/2)(z+l/2)-(n+l/2)(2+l/2)] 

n,n 

_|_ J ^_^-jm e j7r[mgTi + (i/2)(m 2 +g 2 )r2] e 7ri[m(z+2)+g(z-5)] 

m£2Z+l g£2Z me2Z ge2Z+l/ 

(C8) 
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where m = n — n and q = n + n + 1. Applying the Poisson resummation formula 

^exp(-7T(m 2 + bn) = -^J^exp — — (k + b/2ni) 2 (C9) 



a 

k 



for the summation over q in Eq. (1C8I) leads to 

|^(^)| 2 = -^=L ^ e -™ 2 r 2 /2 e ^m( 2 +,-) | ^ eX p[-(7r/2r 2 )(m / -mn-z + ^) 2 ] 

7ri(z— Z+IT2/2) 

+ — e- 7rm2T2/2 e mm{z+ ~ z+Tl) ( ^exp[-(7r/2r 2 )(m / - mn - z + z - ir 2 ) 2 } 



2 



mG2Z 



2 m,m' 

(CIO) 

where (p^' m \z, z) = (n / iT^\m{jz —fz) +m' (z — z)] as in Eq. ( 1C5I) . (Note that the definition 
of the Jacobi theta functions, Eq. (J9j) is also used.) 

Eq. flCSH and Eq. ( 1C10I) tell us that when summing over different winding numbers m 
and m', weighting different winding sector by the sign factor — (— i)( m + 1 )( m '+ 1 ) would result 
in extracting out the a = (1/2, 1/2) sector of the Ising model. It is now possible to calculate 
the boson correlators of Eq. (IC3I) . e.g. 

(d<f)(z)d<f)(z) exp(i0(u>i, w^/2) exp(-i0(w 2 , TS 2 )/2)) a= (i/ 2 ,i/2) 

V}e~ S ^ J2 (-l) (m+1){m ' +1) 4 m ' m ' ) 90 (m ' m ' ) (^)«90 (m ' m ' ) (^) 

m,m' 

x exp(^ m ' m ')(w 1 ,w 1 )/2)exp(-z < /)( m ' m ')(u; 2 ,?Z; 2 )/2) 
-1 

m,m' 

x [d<p i ™' m '\z)d<p i ™' m '\z)(exp(i<p(wi,^ 
+ d(p { Q n,m \z) (d(f){z) exp(i0(to x , tt)i)/2) exp(-20(w 2 , w 2 )/2)) 
+ <94 m,m 5 (2) exp(i0(wx, Wi)/2) exp(-i0(w 2 , w 2 )/2)) 

+ (90(z)90(z)exp(i0(w;i,wi)/2)exp(-i0(w; 2 ,w7 2 )/2))]. (Cll) 

The summation over m and m! can be eliminated by inserting the result of Eq. ( 1C10I) . 



Z n Z 



^E(- 1 ) (m+1)(m ' +1) ^ m ' )ex P^o m ' m ' ) (^/2,^i 2 /2)] 
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together with some necessary differentiations, into Eq. (ICllj) : 



(d(j)(z)d(j)(z) exp(i0(wi, wi)/2) exp(-i(j)(w 2 , w 2 )/2)} a= (i/ 2 ,i/ 2 ) 
\/2t^ 



exp[7r(w 12 - iwia) /8r 2 ][(exp(i0(wi, Wi)/2) exp(-i0(w 2 , w 2 )/2)) 



x {-1^(^/2)1= + I — + ^>™f) |, 1 (.», 12/2 )p 



r 2 



4r| 



H z V 1 {w 12 /2)$ 1 (w 12 /2) 01(1012/2)0! (wia/2)} 

2t 2 2r 2 

- i(8<j)(z) exp(i(f)(wi, wi)/2) exp(— i(f>(w 2 , w 2 )/2)) 
x M^i2/2) U(w 12 /2) + ?r(Wl2 2 ~ Wl2) Mw 12 /2) 

- i(d<fr(z) exp(i<f>(wi, wi)/2) exp(— i(f>(w 2 , w 2 )/2)) 
x $ x {w 12 /2) U(w 12 /2) - ?r(Wl2 2 ~ Wl2) 0i(^i2/2) 
+ (90(z)90(z)exp(#( Wl ,w 1 )/2)exp(-^(«; 2 ,«; 2 )/2))|^ 1 («; 12 /2)| 2 ]. 



(C12) 



Only the correlators of the "free part" of the boson remains to be determined. Its 
propagator is 



ir(z — z) 2 



2t 2 



(<f>(z,z)(f>(Q,0)) = -In 
Note that this propagator is the solution to the modified Green function satisfying 

47T 



-AG{z,z) =A-k5 {2 \z) - 



From Eq. (1013j) . one can obtain 

(exp(i0(wi,Wi)/2) exp(— i<fr(w 2 , w 2 )/2)) = exp 



t 2 



TX[W l2 - W 12 



8r 2 



(C13) 



(C14) 



0'i(O) 
$i(w 12 ) 



1/2 



(C15) 



Since 



{d<j>{z) exp[i(j)(wi,Wi)/2] exp[— i0(w 2 , w 2 )/2]) = exp 



-(0(u;i,wi)0(u; 2 ,w 2 )) 



x (<90(z) exp[i(0(wi, wx) - (f)(w 2 , w 2 ))/2] 



exp 



■(0(Wl,U7l)0(w 2 ,W 2 )) 



x {{d4>{z)^{w u wi)) - (d(/)(z)<f> c (w 2 ,w 2 ))), 

(C16) 
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Eq. (ICT3j) leads to 

(d<fi(z) exp[icf)(w 1 , Wi)/2] exp[-icf)(w 2 , w 2 )/2]) 



2 eXP 



7r(wi2 ~ W 12 f 

8r 2 



gj(g) 

$iOi 2 ) 



1/2 



X 



- Wx) _ $[(z - W 2 ) 7r(Wi2 - Wi 2 ) 

#i(z — wi) di{z — w 2 ) T 2 



(C17) 



Similarly 



(<90O) expf^Ox, Wi)/2] exp[-i0(w 2 , w 2 )/2]) 



7r(wi2 - w 12 f 
8t 2 



^1(^12) 



1/2 



X 



"2 6XP 

^(z -tui) ^(g-wa) | 7r(w 12 - W12) 



r 2 



(C18) 



Lastly 



(d<fr(z)d4>(z) exp(«0(wi, Wi)/2) exp(-i(/)(w 2 , w 2 )/2)) 



exp 



7T 



r 2 



exp 



7r(wi2 ~ W12)' 

8r 2 

7r(Wi2 - ^1 2 ) 2 " 

8r 2 



1/2 



^i(O) 
#1(^12) 
^(0) 1/2 



#l(z - Wi) _ ^[(Z - W 2 ) _ 7r(wi 2 - Wu) 

fi^z — wi) ■di{z — w 2 ) 



T 2 



fii(w 12 ) 



(C19) 



The second term of this equation comes from the fact that from Eq. (1C13I) . (d(p(z)d(j)(z)) is 
actually nonzero. 

Inserting Eq. (jCISl) . 07} . (IUT8|) . and (jUT9l) into Eq. (IUT2|) gives 



{d(f)(z)d(t)(z) exp(i0(wi, tUi)/2) exp(-i0(u;2, w 2 )/2)} a =(i/ 2 ,i/2) 



^1(^12) 



1/2 



4(^2) + ( jj^B) 



giCg - w 2 ) 
- w 2 ) 



(C20) 



Similarly, 



{d(j){z)d(t){z) exp(-i0Oi, wi)/2) exp(i0(w 2 , w 2 )/2)) a =(i/2,i/2) 



^i(O) 
^1(^12) 



1/2 



, ;( - mi2/2 )-I, l( - W2) (gi^-gi^| 



(C21) 
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Eq. ODjl and flC2ll . together with Eq. O and Q leads to the following result stated 
in Sec.III B: 

1/8 



{^(z)a(w 1 )a(w 2 )} a =(i/2,i/2) 



oc 



#1(^12) 



— lt?i) l?]^ — lt? 2 ) 

&i(z — w\) •& 1 {z — w 2 ) 



1/2 



(C22) 
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